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PACS 82 . 35 . Lr - Physical properties of polymers 
PACS 87 . 15 . Vv - Diffusion 

PACS 87. 15. Aa - Theory and modeling; computer simulation 

Abstract. - We study the translocation dynamics of a polymer chain threaded through a nanopore 
by an external force. By means of diverse methods (scaling arguments, fractional calculus and 
Monte Carlo simulation) we show that the relevant dynamic variable, the translocated number 
of segments s(t), displays an anomalous diffusive behavior even in the presence of an external 
force. The anomalous dynamics of the translocation process is governed by the same universal 
exponent a = 2/(2^ + 2 — 71), where v is the Flory exponent and 71 - the surface exponent, 
which was established recently for the case of non-driven polymer chain threading through a 
nanopore. A closed analytic expression for the probability distribution function W(s,t), which 
follows from the relevant fractional Fokker - Planck equation, is derived in terms of the polymer 
chain length N and the applied drag force /. It is found that the average translocation time 
scales as r oc f~ 1 N~~ 1 . Also the corresponding time dependent statistical moments, (s(i)) cx t a 
and (s(i) 2 ) oc t 2a reveal unambiguously the anomalous nature of the translocation dynamics and 
permit direct measurement of a in experiments. These findings are tested and found to be in 
perfect agreement with extensive Monte Carlo (MC) simulations. 
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Introduction. — Recently single molecule experi- 
ments probing single - stranded DNA or RNA transloca- 
tion through a membrane nanopore attracted widespread 
attention [1]. These investigations have been triggered in 
the seminal experimental paper by Kasianowicz et al. [2] 
where an electric field drives single - stranded DNA and 
RNA molecules through the a - hemolysin nanopore so 
that each threading is signaled by the blockage of the ion 
current. By recording the blockage time one can recon- 
struct the whole driven translocation of DNA molecule. 
More recently solid - state nanopores have been used for 
DNA translocation experiment [3,4]. Such pores can be 
tuned in size and are more stable over a wide range of 
voltages, temperature as well as the solvent pH. 

The physical nature of the translocation process is still 
not well understood. The theoretical consideration of the 
translocation dynamics is usually based on the assump- 
tion that the translocation length s (i.e. the translocated 
number of segments at time t) is the only relevant dynamic 
variable which is governed by a conventional Brownian dif- 



fusion process [5-7] . The main predictions for the average 
translocation time r looks as follows. For an unbiased 
translocation t{N) cx a 2 N 2 /D (here a is a polymer Kuhn 
segment length and D is a diffusion coefficient whose N- 
dependence is not well established) whereas the r for the 
driven translocation (when a polymer experiences a chemi- 
cal potential difference A/i between the environments sep- 
arated by the membrane) scaled as t cx Ta 2 N/ (DA/i). 
Here T denotes temperature and we have set the Boltz- 
mann coefficient ks = 1. More recently Kantor & Kar- 
dar [8, 9] have cast doubt on these results by noting that the 
unimpeded motion of a polymer scales as the characteristic 
Rouse time r Rouso cx N 2u+1 where v = 0.588 at d = 3 [10], 
so that TRouse ^ t although the unimpeded motion must 
be faster than that of a constrained chain. Kantor & Kar- 
dar argue [8, 9] that the lower bound of the transloca- 
tion time (which corresponds to an unimpeded motion 
of a polymer) should scale as r(A/i = 0) cx N 2u+1 and 
r(A/i) cx A fl+l/ /A/i for an unbiased and driven translo- 
cation cases correspondingly. They also carried out MC - 
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simulations (using a bond fluctuating model on a 2d lat- 
tice) and the results show that r(A/i = 0) oc A 2 5 and 
t(A/i) oc A 153 /A/i, i.e. at least in the case of the driven 
translocation the theory is inconsistent with MC - simu- 
lation data. A more recent MC and Langevin dynamics 
study [11] reports scaling laws r oc A 15 and r oc A 165 
for relatively short and relatively long chains respectively. 
It has also been questioned [8, 9] whether the transloca- 
tion dynamics is that of normal Brownian motion and 
suggested instead that anomalous diffusion dynamics [12] 
might explain some MC - findings. Nevertheless, there is 
so far no clear knowledge regarding the physical origin of 
such anomalous dynamics. It is also not clear how one can 
make use of the fractional Fokker - Planck equation [12,13] 
which seems to govern this type of dynamics. 

In this paper we come up with a general picture of the 
driven polymer translocation based on our previous con- 
sideration of the unbiased problem [14]. We first sketch 
the mapping of the 3d problem on the Id translocation s - 
coordinate. This leads to an anomalous diffusion in the ex- 
ternal force field which one could quantify in terms of the 
fractional Fokker - Planck equation (FFPE). The solution 
of this equation is then obtained on the interval < s < N 
in a closed analytical form . The subsequent comparison 
of our extensive MC - results with the proper analytical 
expressions shows a very nice quantitative agreement. 
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Fig. 1: How a fold squeezes through a nanopore. The driving 
force / is caused by a chemical potential gradient Afj, = /ii — fj,2. 
(a) The fold of the length s is fragmented into n and s — n - 
parts during its threading, (b) This fragmentation gives rise to 
an effective entropic barrier F(n) with height AE at n = s/2. 



Dynamics in terms of a single translocation coor- 
dinate. — As already noted, the initial 3d problem can 
be rephrased in terms of Id translocation coordinate s, 
and in doing so one arrives at a typical case of anomalous 
diffusion. Recently we suggested [14] that the transloca- 
tion proceeds by successive threading of small fractions of 
the polymer, called folds, which equilibrate fast enough 
compared to the whole chain, and can be considered as 
building blocks of such mapping. In a somewhat different 
context concerning the polymer dynamics the notion of 
folds has been discussed earlier [15]. Figure 1 shows how 
a fold overcomes an entropic barrier caused by a narrow 
pore. 

If the fold is fragmented into n and s — n parts while it 
is threading through the pore then the corresponding free 
energy reads F(n)/T = —sin ft— (71 — 1) ln[n(s— n)], where 
k is the connective constant and 71 is the surface entropic 
exponent [16]. Then the corresponding activation barrier 
which could be associated with the fold threading can be 
calculated as AE(s) = F(s/2) - F(l) = (1 - 7 i)Tms. 

In the force - free case the characteristic time of the fold 
transition from cis - to the trans - side of the membrane 
can be estimated as follows. In the absence of a separat- 
ing membrane this would be the pure Rouse time tR OUS0 oc 
s 2u+1 . The membrane with a nanopore imposes an addi- 
tional entropic activation barrier AE(s) which slows down 
the transition rate. The characteristic time, therefore, 
scales as t(s) = t Rousc (s) exp[A£(s)] oc s 2 ^+ 2 -7i. This 
makes it possible to estimate the mean-squared displace- 



ment of the s - coordinate as 



(s 2 ) oct 2 /( 2 "+ 2 -^. 



(1) 



As a result the mapping on the s coordinate leads to an 
anomalous diffusion law, (s 2 ) oc t a , where a = 2/(2i/ + 
2 — 71). Taking into account that for d — 3 , v = 0.588 and 
71 = 0.680 [17], we obtain a — 0.801. In turn, the average 
translocation time r oc A 2 / Q oc A 2 - 496 . Remarkably in 2d 
where v 2 d = 0.75 and 71 « 0.945 [18], one finds a w 0.783, 
i.e. a is almost dimensionality independent! This explains 
why the measured exponents in both 2d [8] and in 3d [19] 
are so close. The presence of the external force imposed 
on the translocating chain leads to a nonisotropic cis - 
trans - transition of the folds. It can be quantified within 
the FFPE - formalism which was originally suggested by 
Barkai, Metzler and Klafter [20]. 

Fractional Fokker - Planck equation. The for- 
malism of FFPE provides an appropriate technique which 
describes the anomalous diffusion in an external force - 
field. In our case FFPE has the form 
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W(s,t)= Q Dl 



d U'(s) 
ds £ a 



K, 



ds 2 



W(s,t), (2) 



where W(s, t) is the probability distribution function 
(PDF) for having a segment s at time t in the pore, and the 
fractional Ricmann - Liouville operator $D\ ~ a W(s, t) = 
(l/T(a))(d/dt) f*dt'W(s, t')/(t -t'f- a . In Eq. (2) T{a) 
is the Gamma-function, K a is the so called generalized 
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diffusion constant, and £ Q = T/K a is the appropriate fric- 
tion coefficient. In our case the external field U(s) is a 
simple linear function of the translocation s - coordinate, 
namely U (s) — —A[i s, where A/i = \i\ — fi 2 . We consider 
the boundary value problem for FFPE [21] in the inter- 
val < s < N. The boundary conditions correspond 
to the reflecting- ads orbing case, i.e., [U'(s)W(s,t)/T + 
{d/da)W(8,t)]\,=o = and W(s = N,t) = 0. The 
initial distribution is concentrated in so, i.e., W(s,t = 
0) = 5(s — so). The full solution can be represented as 
a sum over all eigenfunctions tp n (s) [22], i.e., W(s,t) = 
exp($(s ) - Hs))En=o T n(t)Ms)Mso) where $(s) = 
U(s)/2T , ip n {s) = exp($(s))<^„(s) and ip n (s) obey 
the equations [(d 2 /ds 2 ) — f (d/ds) + X n .a/K a ] <p n (s) = 
(where / = A/i/T) , and the eigenvalues A„ Q can be read- 
ily found from the foregoing boundary conditions. The 
temporal part T n (t) obeys the equation (d/dt)T n (t) = 
—^n,a qD\ ~ a T n (t). The solution of this equation is given 
by T n (t) = T n (t = 0)E a (-\ niCt t a ) [12] where the Mit- 
tag - Leffler function E a (x) is defined by the series ex- 
pansion E a (x) = J2n=o x k /T(l + ak). At a = 1 it 
turns back into a standard exponential function (nor- 
mal diffusion). Allowing for the boundary conditions 
leads to a transcendental equation for the eigenvalues, i.e. 
-2y/K^/f = tan(^/K^V), where n n = X n , a /K a - / 2 /4. 
This eigenvalue problem has simple solutions in two lim- 
iting cases. For a very weak force fN -C 1 the n n = 
^n,a/K a = (2n + l) 2 7r 2 /47V 2 and the eigenfunctions take 
on the form ip n (s) = y/2/N cos [(2n + \)-ks/2N}. The re- 
sulting solution for W(s, t) at / = reduces to that of the 
force - free case [14]. 

In this paper we focus our attention on the oppo- 
site limit, fN » 1, i.e. when the driving force is 
pretty strong. In this case the eigenvalues spectrum 
reads A„, Q = (/ 2 /4 + n 2 TT 2 /N 2 )K a and the eigenfunctions 
ip n {s) = sin(mrs/iV), so that the resulting solution 

becomes 
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In the limit of strong driving force the translocation 
times are relatively (as compared to the force - free case) 
short and we could use the small argument approximation 
for the Mittag - Leffler function E a (—x), i.e. E a (—x) ~ 
exp[— x/r(l + a)] at i < 1. This makes it possible to 
obtain an explicit analytical expression for W(s, t) which 
can be derived by replacing the summation by an in- 
tegral in eq. (3). In doing so one should use the re- 
lation 2 sin(n7rso/iV) sin(n7rs/A^) = cos[n7r(s — so)/N] — 
cos[n7r(s + sq)/N]. Then one can integrate over n ex- 
plicitly, taking the limit so — ► 0, and finally normalize 
theFPTD: w(s,t) = lim So ^ W(s, t)/ f" W(s,t)ds. This 



yields eventually 



w(s,t) = 



cxp{-(s-/t) 2 /4i} 



Id |erf[/\/|/2] - erf[(/t - N)/2Vt}} 



(4) 



where the dimensionless force / = A/i/T , t = K a t a /T(l + 
a) and erf (x) is the error function. Our further theoretical 
findings are based mainly on eqs. (3) and (4) for PDF. 

First - passage time distribution. — In the chain 
translocation experiment the initial position so can be 
fixed and the distribution of the translocation times is 
actually equivalent to the first - passage time distri- 
bution (FPTD) Q(s ,t) [22]. The relation Q(s a ,t) = 

-(d/dt) f*W(s,t)ds [22] enables to calculate FPTD ex- 
plicitly. Starting from eq. (3) we arrive at the expression 



Q(so,t) = 



N2 t i- a 



OC 

£ 

n=0 



(-D ( - 1) S in(^) 



x E n 



f 2 nV 

1 

4 N 2 



K n t c 
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where the generalized Mittag - Leffler function E a a (x) = 

££=o*Vr(a + *a). 

In the same manner as above we could use the small 
argument approximation for the generalized Mittag - 
Leffler function E a ^ a (— x), i.e. E a . a (— x) ~ (a/T(l + 
a)) exp[— x/r(l + a)] at x <C 1, to obtain the ex- 
plicit analytical expression for Q(so,t). The substi- 
tution of the summation by integration in eq. (5) 
and the use of the relation (— l) n_1 sin(n7rso/-/V) = 
— cos(n7r) sin(n7rso/-/V) = sin[ri7r(l — sq/N)} — sin[n7r(l + 
so/N)] enable finally to obtain for the normalized FPTD, 
lim So ^ Q(so,t)/ J Q(s ,t)dt — > Q(t), the following ex- 
pression 



Q(t) = 



\^l 2 ft 

( 

x exp < — ■ 



r(i + q) 

K a t a 
N - f K - 1 



1/2 



N 2 T(l + a) 



K a t c 



r(l+a) 



4 r(T+S) 



(6) 



As one can see, after normalization the dependence on the 
initial value so — > drops out. It is of interest that FPDT 
given by eq. (6) exactly coincides (at a = 1, i.e., in the 
Brownian dynamics limit) with the corresponding expres- 
sion in the paper by Lubensky & Nelson [23]. It is also 
evident from eq. (6) that the maximum position scales as 
t max cx (N/f)^ a = (N/f) 1 - 25 . Nevertheless, the function 
Q(t) is quite skewed and we will see below that the av- 
erage translocation time r = j tQ(t)dt (which is presum- 
ably measured in an experiment) scales differently. Note 
that Eq. (6) is valid for t < N 2 l a (T(l + a)/2K a ) 1/a , i.e., 
t < 0.47V 2 5 for a = 0.8 which is not a serious limitation 
actually because our translocation times scale as r cx N 1 - 5 
as will be demonstrated below. 
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Statistical moments (s) and (s 2 ) vs. time . 

The recording of statistical moments time dependence, 
(s(i)) = J Q sw(s,t)ds and (s(i) 2 ) = s 2 w(s,t)ds, is 
very instructive (as in the force - free case [14]) for the 
consistency check. Starting from eq. (4) the calculation 
of the first moment yields 



(*(*)) 




exp[ 



-(ft-N) 2 /4t\ 



erf[/Vt/2] - erf[(/t - N)/2y/t\ 



It can easily be shown that in the large time limit (s) 
In the same manner the second moment reads 



< S 2 (t)> = / 2 P + 2i + 2W± 



(7) 



N. 



(8) 



x /f~exp[-/ 2 f~/4] - (/[+ AT) expK/j- W^] 
erf[/Vt/2] - erf[(/t - A r )/2Vt] 

In eqs. (7) and (8) the notations are the same as in eq. 
(4). The detailed check of these relations will be given 
below. Here we only note that for large times 1 < t < N/f 
the exponential terms in eqs. (7) and (8) vanish so that 
to a leading order the moments vary as (s(i)) oc t a and 
(s(i) 2 ) oc t 2a . Again it can be shown that at t — ► oo 
the moments (s(t)) and (s 2 (i)) saturate to plateaus which 
scale like N, and N 2 respectively, as they should. 

Scaling arguments. — The foregoing theoretical con- 
sideration has been based on a rigorous mathematical 
treatment of the FFPE. Before proceeding to the MC - 
check of these findings we put forward some simple scal- 
ing arguments so as to quantify the mean translocation 
time as well as the statistical moments. Let us take the 
average external field energy | (U(s)) | = fN as a natural 
scaling variable. Then the driven translocation rate scales 
as r" 1 = r ~V(/A), where r oc N 2 ^ 2 ^ 1 denotes the 
translocation time in the force - free case [14]. The scaling 
function <f>(x) behaves in the following way: <p(x <C 1) — 1 
and 4>(x » 1) ~ x because at fN > 1 we could expect 
that the translocation rate is proportional to the force /. 
As a result we come to the conclusion that at fN ^> 1 the 
translocation time is scaled as 



r oc y N 2y 



+ 1-71 



(9) 



Taking into account the values for v and 71 given above 
we arrive at the estimations: at d — 3 the translocation 
exponent 9 = 2v + 1 — 71 = 1.496 and at d = 2 the 
exponent 8 = 2v + 1 — 71 = 1.56. This is pretty close to 
the estimation given by Kantor& Kardar [9], 6 = 1.53. 

If we assume that the behavior of (s 2 (t)) (before it hits 
the plateau) follows a power law, (s 2 (i)) oc (ft) 13 , then 



from the correspondence to the scaling law, eq. (9), one 
may estimate (3. Indeed, at the translocation time (fr)^ ~ 
N and the requirement of correspondence with eq. (9) 
yield f3 = 2/(2v + 1 - 7i ). This gives f3 = 1.334 at d = 3. 
In the next section we will demonstrate that this power 
law is in reality too crude. The observed exponents in the 
MC simulation (as well as in the analytical theory given 
above) cross over from a smaller value at very short time 
to a larger one (2a ~ 1.6) at long times in comparison 
with the simple scaling prediction (3 = 1.334 sa 4/3. 

Monte Carlo data vs. theory. — We have carried 
out extensive MC - simulations in order to check the main 
predictions of the foregoing analytical theory. We use a dy- 
namic bead-spring model which has been described before 
[24], therefore we only mention the salient features here. 
Each chain contains N effective monomers (beads), con- 
nected by anharmonic FENE (finitely extensible nonlinear 
elastic) springs, and the nonbonded segments interact by 
a Morse potential. An elementary MC move is performed 
by picking an effective monomer at random and trying to 
displace it from its position to a new one chosen at ran- 
dom. These trial moves are accepted as new configurations 
if they pass the standard Metropolis acceptance test. It is 
well established that such a MC algorithm, based on lo- 
cal moves, realizes Rouse model dynamics for the polymer 
chain. In the course of the simulation we perform suc- 
cessive runs for chain lengths = 16, 32, 64, 128, 256, 512 
whereby a run starts with a configuration with only few 
segments on the trans-side. Each run is stopped, once the 
entire chain moves to the trans-side. Complete retracting 
of the chain back to the cis-side is prohibited by taking 
the head monomer larger than the pore diameter. Dur- 
ing each run we record the translocation time r, and the 
translocation coordinate s(t). Then we average all data 
over typically 10 4 runs. In Fig. 2 we show the PDF Q(t) 
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Fig. 2: First passage time distribution functions at = 128 
and different forces as calculated from MC - data (symbols) 
and the theoretical prediction eq. (6) (solid lines). 
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of a polymer chain with TV = 128 for three different values 
of the drag force, / = 0.5, 0.8 and 1.0. Although the MC 
data is somewhat scattered, especially for / = 0.5, the 
agreement with the analytic expression, eq. (6) is very 
good. Since we set the generalized diffusion coefficient 
K a = 1 and T(l + a) sa 0.931 for a = 0.8, the comparison 
with MC results suggests that a time unit in the FFPE 
corresponds roughly to 500Af CS. 

Using the PDF Q(t), one may determine the MFPT 
(or, translocation times) r which are compared in Fig. 3 
for 16<iV<512 and six values of the drag force /. Ev- 
idently, for both theory and simulation the data collapse 
on master curves fr cx TV 15 , if one scales r with the re- 
spective force, cf. eq.(9). It is seen that the simulation 
data is shifted up by a factor of w 500 which translates 
the MC time into conventional time units. The variation 
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Fig. 3: The average translocation time versus chain lengths. 
The upper line represents the results of MC - simulation, the 
lower refers to the theoretical prediction obtained by the proper 
numerical integration of FPTD eq.(6). Both lines correspond 
to a power - law dependence with exponent 1.5 

of the moments (s), and (s 2 ) is displayed in Fig. 4. Again 
a perfect collapse of the transients is achieved by scaling 
the time with the applied force t — > tf. One can imme- 
diately see that the simple scaling prediction (s 2 ) cx t@ is 
not perfect: for tf < 3 evidently (s) grows with a smaller 
exponent whereas at later times the increase is steeper. As 
mentioned above, this course is very well accounted for by 
eqs. (7), (8). Thus, for tf <C 1 one can readily obtain from 
eq. (7) as a leading term (s) cx W 2 while for 1 < tf < N 
one has (s) cx t a . As indicated in Fig. 4, the observed 
agreement between theory and computer experiment is re- 
markable indeed. Notably this finding suggests that even 
the presence of drag force does not eliminate the anoma- 
lous character of the translocation process as one would 
intuitively expect. This result resolves thus a problem, 
raised initially by Metzler and Klafter [13]. The univer- 
sal exponent a = 2j{2v + 2 — 71) for unbiased threading 
through a pore is not suppressed by the drag force! One 



may thus conclude that the measurement of the number 
of translocated segments with time could provide a means 
for direct observation of anomalous diffusion. 
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Fig. 4: Statistical moments versus reduced time tf from MC 
data and from the analytic results, eqs. (7), (8) for chain length 
16 < N < 256: (a) The first moment (s(t)>: the slope /3/2 = 
2/3 is indicated by a long dashted line, a short dashed line 
denotes (s) oc t a . (b) The second moment ^s 2 (£)^: A long 
dashed line indicates a slope /3 = 4/3, a short dashed line 
denotes (s 2 > cx t 2a . 

Summary. — By solving the fractional Fokker-Planck 
equation for a driven polymer translocation through a nar- 
row pore and deriving a closed analytic expression for the 
probability distribution function W(s,t) to have the seg- 
ment s of the chain in the pore at time t we have demon- 
strated that the translocation process displays all features 
typical for anomalous diffusion. The physical background 
of this behavior is elucidated by scaling considerations. 
The polymer translocation is considered as a squeezing 
of subsequent chain fragments (folds), each being in local 
thermodynamic equilibrium, through a narrow pore. This 
consideration gives rise to an universal scaling exponent 
for anomalous diffusion a — 2/(2v + 2 — 71) so that the 
time t needed for a chain of N segments to move from 
cis to the trans semispace in the absence of drag scales as 
T CX N 2 / a . The presence of external pulling force modi- 
fies this relationship to r cx f~ 1 N 2v + 1 ~n , This principal 
result of the present investigation is unambiguously con- 
firmed by calculation of the mean first passage times (the 
average translocation times) from the derived analytic ex- 
pression for the translocation time distribution function 
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Q(t) as well as by comparison to the results of extensive 
Monte Carlo simulations. We also show that the growth 
of the average number of translocated segments (s) with 
time follows a power law (s(t)) oc t a (for relatively long 
times) which directly displays the anomalous diffusion ex- 
ponent a. Our analytic data also appears to be in perfect 
agreement with the simulation results in a wide range of 
polymer lengths and forces. Thus we have demonstrated 
that the translocation dynamics of a driven polymer chain 
through a narrow pore retains all features of anomalous 
diffusion despite the application of external force. 
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